Unconventional Electromagnetically Induced Transparency in Strongly Correlated 

Quantum Gases 



o 

(N 



C/2 

wo. 



i 

C 

O 

o 



> 
in 
o 

cn 

cn 
o 



X 




H. H. Jen 1 and Daw-Wei Wang 1 ' 2,3 
1 Physics Department, National Tsing Hua University, Hsinchu, Taiwan, R. 0. C. 
2 Frontier Research Center on Fundamental and Applied Sciences of Matter, 
National Tsing Hua University, Hsinchu, Taiwan, R. 0. C. 
3 Physics Division, National Center for Theoretical Sciences, Hsinchu, Taiwan, R. 0. C. 

(Dated: March 15, 2013) 

We develop a general theory to study the electromagnetically induced transparency (EIT) in 
ultracold quantum gases, applicable for both Bose and Fermi gases with arbitrary inter-particle 
interaction strength. We show that, in the weak probe field limit, the EIT spectrum is solely 
determined by the single particle Green's function of the ground state atoms, and reflects interesting 
quantum many-body effects when atoms are virtually coupled to the low-lying Rydberg states. As an 
example, we apply our theory to ID Luttinger liquid, Bose-Mott insulator state, and the superfiuid 
state of two-component Fermi gases, and show how the many-body features can be observed non- 
destructively in the unconventional EIT spectrum. 

PACS numbers: 42.50.Gy, 67.85.-d, 71.10.Pm 



Introduction: In the last decade, systems of ultracold 
atoms open up great opportunities to study many-body 
physics, which may be even not existing in traditional 
condensed matter systems [![. Huge advances in exper- 
imental techniques of tuning interaction strengths and 
generating optical lattice enable studying strongly cor- 
related atomic gases with high controllability. Most of 
the experimental measurements rely on the interference 
pattern of matter wave in different circumstances, for 
example the time-of-fiight experiment, the noise corre- 
lations measurement 2], and the Bragg scattering spec- 
troscopy etc. In-situ imaging is also shown important 
in characterizing the equation of the state and the critical 
properties near the phase transition point 0, H[ . 

Besides of the many-body problem, the quantum con- 
trol and manipulation of light-atom interaction is also an 
extensively studied field in cold atoms. One of the most 
important examples is the Electromagnetically Induced 
Transparency (EIT), which has led to several interest- 
ing subjects such as the dark-state polariton, slow light, 
and induced photon-photon interaction in a cold atomic 
gas (1 0] • EIT experiments in quantum degenerate gases 
H, [9 1 or in Rydberg states [l(| are also recently explored 
by several groups. However, so far the EIT theory is still 
mainly based on the single-particle picture (or within a 
mean field approximation), and therefore cannot clarify 
the non-trivial coupling between light-atom interaction 
and the many-body physics. For example, how is the EIT 
spectrum modified when the ground state is a strongly 
correlated state without single particle excitations? In 
what conditions can a non-trivial many-body effect be 
measured experimentally? Resolving these interesting in- 
terdisciplinary problems should bring significant applica- 
tions toward future study on systems of ultracold atoms. 

In this paper, we develop a general EIT theory in a 



FIG. 1: (Color online) (a) Schematic figure for a standard 
EIT experiment in a strongly interacting atomic gas with the 
A-type scheme: The control and the probe (fii) fields 

couple two hyperfine ground states \g) and \s) with the excited 
state |e) (detunings are A2 and Ai respectively). T is the 
spontaneous decay rate of \e). Filled circles denote atoms with 
mutual interaction U in the ground state, (b) Experimental 
setup for EIT on a 2D array of ID Luttinger liquids in the 
counter-propagating scheme. 



strongly interacting quantum gas, taking into account 
the full kinetics and inter-atom interaction, within the 
weak probe field limit. We explicitly show how the ex- 
perimental EIT spectrum is directly related to the dy- 
namical Green's function, and the quantum many-body 
phenomena is manifested most when atoms are coupled 
to a low-lying Rydberg state due to the sufficient recoil 
energy compared to the decay rate of the excited state 
As an example, we apply our theory to three im- 
portant strongly correlated systems: Luttinger liquid of 
ID Bose gases, Bose-Mott insulator in an optical lattice, 
and the superfiuid state of two-component Fermi gases. 
For Luttinger liquid, we demonstrate a power law de- 
pendence of the EIT spectrum near the resonance, while 
the significant frequency shift and the asymmetric ab- 
sorption spectrum can be identified for the other cases. 
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Our results for the first time demonstrate the strongly 
correlated effects in the EIT experiment, suggesting a 
non-destructive measurement on the quantum gases. 

General EIT theory for quantum gases: We consider 
the conventional EIT setup (A type scheme) as shown 
in Fig. QJa): atoms are initially prepared in the ground 
state \g). The probe (tti) and control (f2 2 ) fields couple 
the ground state to another low energy state \s) and an 
excited state |e) with detunings Ai = o>i — (E e — E g ) and 
A 2 = oj 2 — (E e — E s ) respectively where we set h = 1. 
Here lo 12 are central frequencies of probe and control 
fields, while E g ^ se denote the energies of atoms in differ- 
ent internal states. After using the dipole approximation 
and rotating wave approximation, the total Hamiltonian 
(H tot =H + Hu + Hi(t)) in the rotating frame becomes: 



agonalized to be H se — ^ k e_(k)d k dk + e+( i Kjo k Ok 
Here the new eigenstates are: dk = cos ^kek+ki + 
sin^kSk+kj-ka and b k = sjnj^ek+ki - cos^kSk+kx-ka , 
where cos k = y/[e+(k) - eo.k+ki + Ai]/[e+(k) - e__(k)]. 
The associated eigenvalues are: e±(k)=— Aj + [A 2 + 
£o,k+ kl ± V(A 2 -eo,k+k 1 ) 2 +4n|]/2, where A 2 = A 2 + 
£o,k+k r , £o,k = k 2 /(2m) — (j,, and k r = ki — k 2 as the re- 
coil momentum. A phenomenonological spontaneous de- 
cay rate (T) of the excited state can be added by re- 
placing eo,k+ki with eo,k+ki — iT. We also assume neg- 
ligible dephasing rate between the two hyperfine ground 
states. Note that when considering standard D2 transi- 
tion, T s» 6 MHz, is much larger than the average atomic 
kinetic energy, and therefore the quantum many-body 
effect may not be easily observed. 

Defining the polarization operator, P(r, t) 



H = H(^7-M)c k Ck + 5Z[(A 2 - Ai)J k J k do$t(r,t)i> g (T,t) + h.c], where #, e (r,t) is the 
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Here 6k and/or dk denote the atomic field operators in 
the momentum space. Hq includes all the single particle 
part without the probe field for atomic mass m, chemical 
potential \i, and the Rabi frequency of plane wave con- 
trol field fi 2 . All the inter-atom interactions are included 
in Hjj with the short-range interaction U c d (extension to 
a certain long-ranged interaction is straightforward) . Fi- 
nally, Hi denotes the effect of the probe field with VL\^{t) 
being the slow-varying Rabi frequency of wavevector k. 
V is the quantization volume and ki, 2 are momenta of 
the probe and control fields respectively. 

Different from the standard EIT theory, here we keep 
the full interaction effects but treat the probe field part, 
Hi, as a time-dependent perturbation. In other words, 
in the leading order limit of a weak probe, all atoms are 
in the ground state (\g)) and none in states \e) and |s). 
As a result, the interaction between \g) and the other 
two states are nothing but a background energy shift in 
the single particle energy, i.e. E s / e — > E s / e + nU g:S / e , 
where n is the total particle density. The mutual inter- 
action between \s) and |e) are the second order effect and 
hence negligible. Similar treatment can be also easily ap- 
plied to long-ranged interaction if Rydberg states of high 
principal quantum number are involved jl 31 ] . 

Therefore, H = Hq + Hjj can be separated into 
two decoupled parts: H = H 9 + H se , where H 9 in- 
cludes all the ground state kinetic energy and inter- 
action energy terms, while H se includes all the sin- 
gle particle terms of states \s) and |e), and can be di- 
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field operator, and do is the dipole moment, we calculate 
its variation with respect to the probe field (-Hi) via 
the linear response theory [1]: S(P(r,t)) = i J_ dt' 

H{^G\[Hi, H (t'),P H (r,t)]\^ G )H, where Hi,a(t) and 
Pn(j,t) are corresponding operators in Heisenberg 
picture, and I'I'g).^ is the ground state of the unper- 
turbed Hamiltonian H = H + H v = H 9 + H se . In the 
momentum-frequency space, the electric susceptibility 
can be obtained to be (see Supplemental material poj): 

do 
V 
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Here G < (k, ui) is the Fourier transform of the Green's 
function at zero temperature [f|: G < (0, 0;r, t) = 
Ti H ( x i'G\i'l(r,t)il> g (0,0)\^ G )HO(t), where the upper 
(lower) sign refers to bosons (fermions). We note that 
the above result is very general in the limit of weak 
probe field, and all the many-body effects are hidden in 
G < (k, oj). However, as pointed out in Ref. 21!], such 
many-body effect can be identified more easily in the 
counter-propagating (rather than the co-propagating) 
setup of the EIT experiment, since the recoil momen- 
tum of the former is manifested more significantly in the 
quantum statistics. In the rest of this paper, therefore, 
we just consider the counter-propagating scheme in order 
to measure the many-body effect from the EIT spectrum. 

EIT in a weakly interacting condensate: We first con- 
sider a weakly interacting condensate as an exam- 
ple. By separating the fluctuation from the conden- 
sate: gk{t) = \/iVc<$k=o + 8gu.(t) with N c being 
the condensate particle number, we apply Bogoliubov 
transformation to calculate the single particle Green's 
function. The final EIT spectrum can be separated 
into the condensate (C) and the non-condensate (NC) 
parts: x(q,w) = xcM + Xnc(<1,u), where xcM = 
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, vl [cos 2 



d "c[cos 2 0o /(w + £-(0)) + sin 2 o /(a; + e+(0))], and 
Xnc{u) = doV- 1 £ k#0 vl\ 
sin 2 k /(w + e+(k) + e(k))]. 
sate density, e(k) = (e 2 (k) 



e_(k) + e(k)) 
N c /V is conden- 
e 2 .) 1 / 2 is the phonon exci- 
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tation energy, and u 2 . = sinh 2 6*(k), where tanh2#(k) 
£2/ei(k), ei(k) = k 2 /2m + n c U gg , and e 2 = n c U gg . Note 
that the non-condensate part is contributed by all mo- 
mentum channels due to the quantum depletion, consis- 
tent with results from the dark state approach (2l| . 

EIT in Luttinger liquid: Now we study EIT in a Lut- 
tinger liquid (LL) as shown in Fig. QJb) , which is a very 
general ID effective model @, [22[ and has no conden- 
sate even at zero temperature. Since all the elementary 
excitations of a LL are collective, it is therefore inter- 
esting to investigate how the probe field propagates in- 
side such strongly correlated system. The single particle 
Green's function can be exactly calculated by bosoniza- 
tion method [2[ (see Supplemental material [20(): 



iG< L {x,t) 



1/2k 



[x 2 + (a + ivt) 2 ] 1/4k ' 



(3) 



where n is the Luttinger parameter, v is the phonon ve- 
locity, and a -1 is the system-dependent momentum cut- 
off. 1 < k < 00 for a short-ranged repulsive interaction, 
while k can be smaller than one if the interaction is long- 
ranged. 

Before inserting G^ L (x,t) into Eq. @ for numerical 
calculation, it is more instructive to first consider ana- 
lytical results in the strong control field limit. We have 
(see Supplemental Material [20|): 



X LL(oj)=d n[-) r(l-— ) 
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where cos < 



V(e++Ai)/(e+-e_), and e± = -Ax 



(A 2 ± y / Af+4fi|)/2. r(x) is the gamma function. Note 
that in this limit, there is no momentum (q) dependence, 
and a non-trivial power law dependence appears at the 
resonance frequency due to the LL nature. This power 
law dependence also suggests that the zero of the dis- 
persion relation does not correspond to the transparency 
point (see Supplemental material 20] ) , which reflects one 
of the most important features in the EIT spectrum of 
LL. 

In Figs. Efa) and (b), we show the full numerically 
integrated EIT absorption and dispersion profiles (from 
Eq. © and (0) of a ID 87 Rb atoms (see the caption 
for parameters) in the counter-propagating scheme. It is 
easy to see that the power law dependence of \ll makes 
the dispersion highly asymmetric about the resonance 
points and the absorption depth (a) becomes larger at 
the transparency point (A*) in the strong interacting 
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FIG. 2: (Color online) The EIT profiles for a Luttinger liquid 
of 87 Rb atoms in the counter-propagating excitation scheme. 
We take a static probe field (q, w = 0) and a resonant control 
field (A2 = 0) with Rabi frequency Q.2 = 5 T. The excited 
state is chosen as low- lying Rydberg transition of I24P3/2) 
with the spontaneous decay rate, T = 28.3 /j,s. The effec- 
tive ID interaction strength can be tuned by the confinement 
resonance of ID tubes [2J]. (a) Absorption (Re[ix]) and (b) 
dispersion (Im[ix]) profiles for n — 1 (solid) and 10 (dash). 
The horizontal line guides the eye to the zero in (b). (c) shows 
the group velocity v g (□), and the absorption depth a (0) at 
the transparency point (A*) as a function of k. For con- 
venience, we take the phonon velocity «=4.3 mm/s and the 
cut-off a=0.12 fiia respectively within a typical experimental 
regime [25( |. 



regime (smaller k), in contrast to the standard EIT pro- 
files in the weak interacting limit (rt 3> 1). This is be- 
cause when a ground state atom is excited by the probe 
field, the transition matrix elements are composite of all 
the collective excitations of different energies, weighted 
by the density of states in a power-law distribution. In 
Fig. HKc), we show the group velocity of light (propor- 
tional to the slope of dispersion at the transparent point, 
where the absorption is the smallest) as a function of k. 
One can see that for a stronger interaction (smaller k), 
the group velocity becomes smaller due to steeper dis- 
persion relations, but the absorption is also larger, mak- 
ing the slow light propagate less effectively in a LL. On 
the other hand, such unconventional EIT spectrum also 
implies a sensitive measurement about the many-body 
properties. 

EIT in a Mott insulator: Now we consider the Mott- 
Insulator (MI) of strongly interacting bosons in a 3D 
optical lattice, which can be well described by a single 
band Hubbard model (HM) j^. The ground state field 
operator is expressed as $ 9 (r, t) = ^R5R.(^)u ; R,( r ) where 
%(r) is the Wannier function, and gn(t) is the field op- 
erator at site R. When deep inside the MI state, we 
can use the three-state model d, 0] to control the small 
number fluctuation above the Fock state, and obtain the 
Green's function at zero temperature (see Supplementary 
Material 0): iG<(k,t) = n Q £ R |i2} R (k)| 2 e- ieh ( k )'0(i), 
where uq is the integer filling fraction, WR(k) is the 
Fourier transform of wnfr), and e^(k) = eo(k)/2 + 5/i + 
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FIG. 3: (Color online) EIT profiles of a unit-filling Mott- 
insulator with 87 Rb atoms loaded in a 3D optical lattice (lat- 
tice constant d — 426 nm). The corresponding parameters 
(U, J) in unit of the recoil energy Er are (0.5, 0.007) and 
(1.05, 10~ 4 ) for Vo/E R = 15 (solid-blue and dash-red) and 40 
(0 and □) respectively. Absorption (Re[ix]> solid- blue, 0) and 
dispersion (Im[ix], dash-red, □) profiles are plotted for vari- 
ous control field strengths (a) Q2 = 2, and (b) 8 V. All other 
parameters are the same as Fig. 2. 



w(k) is the hole excitation with w(k) ps U/2 when deep 
inside the Mott state. Here eo(k) = 2 J^ a=1 cos(k Q <i), 
and Sfi = —3 J for a 3D square lattice of lattice constant 
d. J and U are the tunnelling amplitude and on-site in- 
teraction. As a result, the susceptibility from Eq. © 
becomes 



\ \//(q.~'i = ^-^2\wn(k)\ 2 

k 

sin 2 </> k+q 



V 



COS" k 4 



+ 



uj + e h (k) + e+(k + q) 



uj + e h (k) + e_(k + q) 
(5) 



Compared to the standard EIT spectrum, we find that 
Xmi is contributed by EIT in different momentum states, 
weighting by the Wannier function of the optical lattice. 
When the momentum distribution is large in a strong 
lattice strength, the contribution from different momenta 
makes the EIT spectrum much broadened. 

In Fig. |3l we demonstrate the EIT spectrum inside 
the MI state of 87 Rb with unit filling for different con- 
trol fields and lattice strengths. As discussed above, the 
momentum distribution of Wannier function lead to an 
asymmetric and inhomogeneous broadening of the ab- 
sorption profiles near the resonance. It makes the EIT 
window less transparent when the control field is weak 
(Fig. St &)), while it becomes less effective in the large 
field limit. In other words, by measuring the absorption 
depth one may also determine the momentum distribu- 
tion of the underlying many-body system. 

EIT in a BCS superfluid state: Finally we investigate 
the EIT in a BCS superfluid of two-component Fermi 
gases, using similar setup [12] by coupling one of the 
atomic ground state to a low-lying Rydberg state 
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Without losing generality, here we just presume the exis- 
tence of a superconducting gap, Abcs, not specifying the 
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FIG. 4: (Color online) EIT profiles of a BCS superfluid 
state of two-component Fermi gases ( 40 K). The excited state 
is chosen as low-lying Rydberg transition of |21P 3 / 2 ) with 
r _1 = 25.4 /is. We fix the chemical potential /1 = cf = 6.7 
hF with atomic density n ~ 10 14 cm -3 and resonant control 
field (A2 = 0). Absorption (Re[i%], solid-blue, 0) and dis- 
persion (Im[ix], dash-red, □) profiles are plotted for pairing 
gap As/cf = 0.1 (solid-blue and dash-red) and 1 (0 and □), 
while the control field strengths are (a) Q.2 = 2 and (b) 10 Y. 
(c) shows the group velocity v g (+) and the absorption depth 
a (O) at the transparency point as a function of Arcs- (d) 
shows the transparency positions Aj (x). 



pairing mechanism. One of the two components (pseudo- 
spin up) is our ground state (\g)) in the EIT experiment, 
and the other is not involved in the EIT process. The 
Green's function can be calculated by transforming to- 
ward the Bogoliubov quasi-particles (a and 0): <7k,t = 
cos^kAk + sin#k/3l k , where sin 2 6*k = (1 — £k/-Ek)/2. 
Here = \J &? BCS + £k i s the excitation energy with 
£fc = k 2 /(2m) — /1 and /i being the chemical potential. 



We find that 
plementary Material [2 



sin^ 6» 2 e- i£k '6»(t) (See Sup- 



and then 



Xscs(q,w) = y ^sin 2 6» k ( 



cos 2 </> k +q 



£ k + £-(k + q) 



+ 



£ k + e + (k + q) 



(6) 



Note that in the limit of zero pairing gap, Abcs — 
0, sin 2 #k becomes a step function at Fermi energy, and 
hen ce x rcs becomes the result of a noninteracting Fermi 
gas 0. 

In Figs. @|a) an d (b), we plot the EIT profiles of the 
BCS superfluid phase for 40 K in the counter-propagating 
excitation schemes. Different control field strengths and 
different gap amplitudes are shown together for compar- 
ison. One can see that when the control field is turned 
smaller, the recoil energy becomes more important so 
that the presence of Fermi sea makes the EIT spectrum 
disappeared, similar to the noninteracting fermions [2lj |. 
When the control field becomes stronger, the quantum 
degeneracy becomes less important so that the standard 
EIT features appear, while an asymmetric broadening 
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appears when the gap is opened. Interestingly, we find 
the group velocity is larger for a larger Abcs (see Fig. 
4(c)), different from results of a LL. In (d) we show that 
the transparency points A* scales as the gap energy, and 
therefore can be used to measure the pairing gap. 

In summary, we develop a very general and analytic 
theory to investigate the EIT spectrum of quantum de- 
generate ultracold atoms. We show that the spectrum is 
solely determined by the single particle Green's function, 
and has most significant quantum many-body physics 
when the atoms are coupled to a low-lying Rydberg state. 
In the regime of strong interaction, the absorption pro- 
file have a highly asymmetric inhomogeneous broadening 
with a frequency shift in detuning. Our results suggest a 
new non-destructive method to investigate the strongly 
correlated physics from the well-developed EIT experi- 
ment. 

We appreciate fruitful discussions with Ite A. Yu, Y.- 
C. Chen, M. Fleischhauer, and G. Juzeliunas. This work 
is supported by NSC and NCTS grants in Taiwan. 
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Supplemental Material: 
Linear response theory for the EIT spectrum 

We apply the linear response theory [lj to the electric 
susceptibility of EIT. For the perturbation Hamiltonian 
Hi(t), the linear response of any operator O under the 
ground state \^g)h of the unperturbed Hamiltonian H 
is {%= 1) 



S(0(t))=il dt' H ^ G H hH (t'),0 H (t) 



H 

(S7) 



where Hi t n{t) — e H\(t)e~ is the operator defined 
in the Heisenberg picture (with subscript H), and the 
same as for Oji(t). 

To calculate the electric susceptibility of EIT, we con- 
sider the variation of the polarization operator, 



P(k,t) = / drP(r,t)e- jk r 

= d J2 [^-k+qW&lW + 5 q We k+q W] , (S8) 
q 



and put the above into Eq. (IS7[) . we have 



<5(P(q,m =i / dt'9(t-t') ( * c 



k.k' 



c. e 



-iHt' 



d e i6t (e + _ q+q , (t)Sq' («) + ffi (t) 



J k,k',q' 



-ifft 



^'eU(t')^(i')e- tHt ', e ,ffi ^(t)e, +q '(t)e- 



(S9) 



where in the last line we have neglected the terms of 
ge^ g and g^eg^e, which are vanishing in the expectation 

value of ground state wavefunction, \^g)h- 

For the unperturbed Hamiltonian (or probe-free 



Hamiltonian, H — Ho + Hjj), all atoms are in the atomic 
ground state \g), and states |e) and \s) are empty. There- 
fore, further expansion of the above equation gives 




Since the excited state |e) and the second ground 
state \s) are coupled by the control field, we 
then replace e k -|_ kl (i) with the new eigen-bases 
and 6k with corresponding eigenenergy e±(k) respec- 
tively. Here a k = cos0 k ek+ki + sin^kSk+ki-ks 
and 6 k = sin k e k +ki - cos0 k Sk+ki-k 2i where 
cos0k-v/[ £ +( k ) - £ o.k+ kl + Ai]/[e+(k) - e_(k)]. The 



associated eigenvalues are: e±(k)=— Ai + [A2 + £o,k+ki ± 
V / (A 2 -eo,k+k 1 ) 2 + 4^2]/2, where A 2 = A 2 + e 0) k+k r , 
eo.k — k 2 /(2m) — /1, and k r = ki — k 2 as the recoil 
momentum. A phenomenonological spontaneous decay 
rate (r) of the excited state can be added by replacing 
£o,k+ki with eo,k+ki — *r. After shifting the momentum, 
q — > q ki and k — > k — ki , we have 
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idp 
V 



dt>d(t-t>) ]T {n lik (f) H <*GUt(t) 

1 k,k',q' 

cos0 k+k <a k+k ,(f') - sin0k+k'& k+k /(*')l &'(*') I*g) h - ft.c(q -> -q)| 
= ^ j°° dt'9(t - f) {fii,,(f )^k(* - [cos 2 «^ k+q e---( k + t >)(*- t ' 

./-oo k 

- ft,.c.(q -> -q)j, 

r 



cos</> q+q <a q+q '(f) - sin0 q+q '& q+q '(t) 



sin 2 ^ k+q e- l£ + (k+q)( ^ t ' ) 



(Sll) 



where we have defined the ground state correlation 
function as F k (f -t')5^= H {^a\gt{t)gk(t')\^G)n- Use 
the Fourier transform 6(P(q,u}))= dtS(P(q,t))e~ iuJt 
, and define the electric susceptibility as x(q> w ) = 
S(P(q,uj)) /fli^q(uj) where f2i iq (w) is the Fourier trans- 
form of Qi !q (t), we have (after the change of variables, 
t = f + t",'dt = dt") 

x(q,w) 



= 1fE/ dtd(t)F k (t) 



V 



+ sin 2 </> k+q e 



cos 2 </»k +q e- l ^+ e -( k+t i» f 



-i(w+e + (k+q))t 



(S12) 



Finally, since the correlation function F k (f — f') 
is related to the time-ordered retarded Green's func- 
tion via 9(t)F k (t) = J dre ikr iG < (0,0;r,t), where 
iG<(0,0;r,t) = ff (* G |y}t(r, f)^(0, 0)|* G ) ff 0(f), 
rewrite above results to be 



we can 



id sr-^ 

x(q.w) = yL 

k 



'k+q 



dwiG<(k,w) 



'k+q 



w — uj — e_ (k + q) uj — u) — e+ (k + q) 



(S13) 



where G < (k, w) is the Fourier transforms of the Green's 
function iG < (0, 0; r, f). The eigenvalues e±(k) and <^ k 
have been defined before. 

It is instructive to simplify the above further in the 
strong control field limit, i.e. O2 is much larger than 
atomic kinetic and interaction energies. It is straightfor- 
ward to have following leading order results: 



X(w) 



do \ - 
V 4^ 



duiiG < (\i., uj] 



cos" 



uj — uj — e_ 



Single particle Green's function of a Luttinger liquid 

Here we demonstrate how to derive the single particle 
Green's function [l[ for a Luttinger liquid (LL). First, 
we use the density-phase representation of low-energy 
bosonic field operators, ty g {x,t) — yjp(x, t)e~ l ^ x,t \ for 
the Luttinger liquid [2J, where p{x, t) and (f)(x, f) are den- 
sity and phase fluctuation operators. The single particle 
Green's function becomes 

iG$ L {Q,Q;x,t) 

~ s/p(x)p(0) / e -ikx,t) e i$(0,0)\ e ^ 



-i<t>(x,t) A<t>(0,0) 



e(t), 



(S15) 



where we have used the fact that the density fluctuation 
is suppressed and negligible for a repulsively interacting 
ID gas [2j, and therefore only phase fluctuation is kept for 
the low energy effective theory. 0(f) is a Heaviside step 
function. The above equation can be simplified by using 
e A e B = e A + B e l A > B }/ 2 if A and B are linear combination 
of bosonic operators and [A, B] is a complex number. 
We can also apply (e A ) = eS A ^ 2 for the expectation on 
a bilinear Hamiltonian of A. In the low energy limit, 
it has been shown that the effective Hamiltonian of a 
ID bosonic gas can be described by a LL model, where 
the phase operator (f>(x, f) can be calculated within the 
periodic boundary condition [2J, 



4>{x) 



2n 



qLn 



1/2 



2 1-? 

x e - Q l 9 l/ 2 sgn(g) \e lqx b(q) + e- iqx fr(qj\ , (S16) 



sin 



uj — uj — e + 



(S14) 



where cos0 = y/ (e+ + Ai)/(e + — e_), and e± = — Ai + 
(A 2 ± ^/A| + 4Qf) /2. Note that xM obtained in this 
limit has no momentum (q) dependence. 



where b(q) is the bosonic eigenstate field operator of the 
effective LL model. A positive length scale, a, is intro- 
duced as a cutoff length scale for the convergence of in- 
tegrals, and L is the system size. Here the LL parameter 
is denoted as k, and the dispersion is oj(q) = \q\v with v 
being the phonon velocity. Note that the density-phase 
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representation in a LL is valid in the long wavelength 
limit q <C p ■ The exact values of k and v should be de- 
termined by a more microscopic calculation or from the 
experimental measurement. 

As a result, the single particle Green's function of the 
LL model can be calculated to be 



iG£ L (0,0;r,t) 

~ e [fa,t),k0)]/2 nexp 
^ e [0(*,t),0(O)]/2 nexp 



T 



(x,t)-<j>(0) 8{t) 



1 x 2 + (a + ivt) 2 ' ( 
~Tn ln a 2 > 9 ®> 



(S17) 

where T is the time-ordered operator. This correlation 
function from LL model describes the behavior of quasi- 
long range order that has an exponent proportional to 
the interaction strength. The logarithmic function inside 
the exponential function is derived from (use uj q = \q\v 
andEq. (IST61) ) 



T 
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E 



2tt 




2 — 2e l ( qx ~\ q \ vt ^ 


qLn 







aq 



2k J q q 

6(t) , x 2 + (a + ivt) 2 
In - . 



2 _ g^O-ft) _ e -iq(x+vt) 



(S18) 

2k a z 

which can be also found in Ref. The commutation 
relation in the prefactor of the Green's function in Eq. 
(|S17|) can be also calculated to be 



<t>(x,t),<l>(0) 



1 

4K 



-a\q\ 



J(qx-\q\vt) 



-i(qx-\q\vt) 



1 j [a + i (x — vt)] [a — i(x + vt)} 
~ AK ° S [a-i(x- vt)} [a + i (x + vt)} ~ ' 



(S19) 



where in the last line we have taken the limit a <g. \x\,v\t\. 

Finally, the full dynamical correlation function of a LL 
model (which is in general valid only in the low-energy 
and long wavelength limit of a ID Bose gas) is then de- 
rived as 



iG< L (0,0;x,t) 



,l/(2«) 



[x 2 + (a + ivt) 2 } 1 /^)' 



(S20) 



The electric susceptibility of Luttinger liquid 

Here we proceed to calculate the electric susceptibility 
(i.e. EIT spectrum) of a LL. It is more instructive to 



consider the case of large control field limit (i.e. Eq. 
(|S14p ) so that by inserting the result of Eq. (|S20|) . we 
have 



Xll{u) 

_ ido 
~ 2ir 



dk / dt 

-oo Jq 



dx- 



ne ikx a l/(2 K ,) 



x 2 + {a + ivt) 2 } 1 /( iK ) 



cos 2 6e 



idon 



-i( W +e_)t + sm 2 ^ e -i(a;+ e+ )i 
a V(2«) 



dt 



o (a + iwi) 1 /(2K) 



cos 2 Se- 1 ^-^ 



-sin 2 6e- i{uj+t + )t 



= *„(-) r ( i--)x 



cos 



sin 



(S21) 



(a; + e_) 1 - 1 /(2«) ( w + e+ )l-l/(2«)_ 
where in the last line we have used the integral property 



dt- 



o (a + ivt) 1 



5 a "/ v (iv)- b (iGj)- 1+b T(l-b,—), 

V 

(S22) 



under the conditions that Re [a] > 0, Im[a;] < 0, and 
Re[u] > 0. r(s,x) is the incomplete gamma function 
which becomes the gamma function T(s) when a — > 0. 

Eq. (|S21[) indicates a nontrivial power law dependence 
in the EIT spectrum, as shown in Fig. 1 from strong to 
weak interacting regimes (k = 0.6-10)0. The standard 
(noninteracting) EIT spectrum is similar to the weak in- 
teracting one (Fig. 1(d)) where the zero of the dispersion 
relation coincides with the transparency point (Ai = 0), 
and symmetric (anti-symmetric) absorption (dispersion) 
profile is retrieved around Ai — 0. When the interaction 
becomes stronger (smaller k), the zero of dispersion re- 
lation moves away from the transparency point, and the 
EIT is destroyed for even stronger interaction. Interest- 
ingly, the EIT profile for k = 1 case (hard core boson 
limit) shows an inversion symmetry between absorption 
and dispersion profiles around Ai = 0, very different 
from the standard EIT spectrum in the noninteracting 
limit (k = 10). 



Single particle Green's function for Mott state 

Here we show how to derive the single particle Green's 
function for a Mott state. At zero temperature, and 
in the deep Mott-insulating limit of average n particle 
per site, we may assume only small number fluctuations 
about no so that the Hilbert space can be truncated to 
three particle numbers per site only: nn — 1, no, and 
no + 1. Within such three-state model 0,13], the original 
bosonic field operator at site R can be re-expressed to 
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FIG. 5: (Color online) EIT profiles for a LL of 87 Rb atoms in 
the large control field limit from Eq. (|S21[) . We take a static 
probe field (q, uj = 0) and a resonant control field (A2 = 0) 
with Rabi frequency Q2 = 5 F. The excited state is chosen as 
low-lying Rydberg transition of I24P3/2} with F _1 = 28.3 fj,s. 
The absorption (Re[ix], solid- blue) and dispersion (ImJixL 
dash-red) profiles for (a) k — 0.6, (b) 1, (c) 2, and (d) 10. 
The horizontal line guides the eye to the zero. 



be: g R (t) = Vno+T£[ )R (t)io,H.(t) + ^tl >R {t)i- 1<n (t), 
where i±i,R and io jR are the raising and lowering op- 
erators. The conservation of total number of particles 
provides an additional constraint on the Hilbert space: 

E q=± i,o^,r(*)*«,rW = 1 - 



Since we are interested in the deep Mott insulator 
regime with little number fluctuation, we can apply the 
conservation of particle number shown above to have 
the following approximation ((J Q: £o,B,(i)> *o R (£) — 

1 - |*1 r (*)*i,r(*) - |*Li r(*)*-i,r(*)- As a result, the 
single particle correlation function can then be easily ex- 
pressed to the quadratic order of small fluctuations to 
be 



<2rW.9R'(0)) =(no + l)<tl lR (t)&,R(*)4,R'( )*i.R'(0)> + n Q (il R (t)L hR (t)il 1RI (0)i , R >(0)) 

+ y/no(n + 1) [<*I r(*#>,r(*)£i R'(°)*o,R'(0)) + r(*)*-i,r(*)*1 r'(0)*i,r'(0)) 



<no + l)^l, R («)ti,R'(0)) + no(*-i,R(*)*Li,R'(0)) + VM^TT) 
x (*lH(*)*li.R'(0)) + <*-i,R(*)*i,a'(0)) ■ 



(S23) 



Within the MI state, it has been shown that we can di- 
agonalize the effective Hamiltonian [6j in the three state 
model, and change the lowering and raising operators 
by the quasi-hole and quasi-particle excitations: j3h and 
/3 P , i.e. Li, k (f) = -B(k))8j jk (t) - A(k)&, k (t), and 

*l,-k(t) = ^0<)$,k(*) + B(k)^,k(t), where 



A(k)=cosh(^), B(k) = sinh(^), 



tanh(Z?k) = 



-2e (k)^/n (n + 1) 
f7-e (k)(2n + l) : 



(S24) 



e (k) = 2J cos(fc a a). 



The corresponding eigenenergies of the particle and hole 
excitations are e P! / l (k) = =F[eo(k)/2 + Sfx] +cD(k), where 
w(k) = ^C/ 2 - J7e (k)(4nn + 2) + e§(k)/2, and <J/x = 
— 3 J for a 3D cubic lattice jfjj. 

As a result, we can easily calculate a correlation func- 



tion as following: 

(*U(*)*1,R'(0)> 



(S25) 



Similarly the single particle Green's function in Eq. (|S23j) 
becomes 

2 



(0)) =^3 E [v/R + IW) - 0*A(k)' 



x e -ik-(R-R') e -ie h (k)t_ 



(S26) 



In the deep Mott regime, we have [/ 3> J, and hence 
A(k) -> 1, and (ff^(t)sjy(0)) = n e'^ k ^ h S R ^. The 
single particle Green's function for the Mott state could 
be further simplified to be (use V> s (k) as the Fourier 
transform of tp g (r) = Xi R 5RWR(r)) 

iG<(k,t) = (^(k,*)^(k,0)Wi), 



= £ ^(k)d)R,(k) (ak(t)flR/(0)) 9{t) 

R,R' 

= ]T|£ R (k)| 2 n e-^ {k)t 9(t). (S27) 



R 
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Note that only hole excitation energy appears because 
the Green's function we need for EIT spectrum is time- 
ordered, i.e. an atom is excited from the ground state 
(\g)) to the excited state (|e)) by the probe field, leaving 
a hole excitation inside the strongly interacting system. 

Single particle Green's function for superfluid case 
of two-component Fermi gases 

Here we show how to derive the single particle Green's 
function for a BCS superfluid state of two-component 
Fermi gases. The single particle Green's function in this 
context is defined as 

~*G< cs (k, t) = (gl t (t)g Kt (0)) 6(t), (S28) 

where we define the original ground state (in the EIT 
A scheme) to be spin up. The other state of spin down 
is assumed not directly involved in the EIT experiment. 
To evaluate the above expectation, we can express the 
ground state in terms of Bogoliubov quasi-particles (de- 
noted by a and (3) of Cooper pairs in the superconducting 
state, 

5k,t = cos #kQ!k + sin k/ SL k , (S29) 

where sin 2 6 k = (1 - &/£ k )/2. Here E k = ^A 2 + 
is the excitation energy of the quasi-particles with = 



k 2 /(2m) — [i. /i is the chemical potential and is deter- 
mined by the particle density. The Green's function can 
then be obtained to be 



-iG< cs {k,t) = ^cos6 k ale iEkt + sme k /3^e- iEkt ^ x 
(cos0 k a k + sin0 k /3l k )^0(i), 
= sin 2 6» k e- l£;k *6>(t). (S30) 
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